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ABSTRACT 

Planck results have revealed that the electric dipole emission from polycyclic aromatic hydrocarbons (PAHs) is the 
most likely explanation for anomalous microwave emission that interferes with cosmic microwave background (CMB) 
radiation experiments. The emerging question is to what extent this emission component contaminates to the polarized 
CMB radiation. We present constraints on polarized dust emission for the model of grain size distribution and grain 
alignment that best fits to observed extinction and polarization data. Two stars with a prominent polarization excess 
at A = 2175 A, HD 197770 and HD 147933-4, arc chosen for our study. For HD 197770, we find that the model 
with aligned silicate grains plus weakly aligned PAHs can reproduce the 2175 A polarization feature; whereas, for HD 
147933-4, we find that the alignment by silicate grains only can account for that feature. The alignment function of 
PAHs for the best fit model to the HD 197770 data is employed to constrain polarized spinning dust emission. We find 
that the degree of polarization of spinning dust emission is about 3.2% at frequency v = b GHz and rapidly declines 
to below 1.8% for v > 20 GHz. We also predict the degree of polarization of thermal dust emission at 353 GHz to be 
Pem = 9-9% and 14.5% for the HD 197770 and HD 147933-4 stars, respectively. 

Subject headings: cosmic microwave background-diffuse radiation-dust, extinction-radiation 
mechanisms : non- thermal 



I 1. INTRODUCTION 

Cosmic Microwave Background (CMB) experiments 

] (see Bouchet et al. 1999; Tegmark et al. 2000; Efstathiou 
2003; Bennett et al. 2003) are of great importance for 
studying the early universe and its subsequent expan- 

[ sion. Precision cosmology with Wilkinson Microwave 
Anisotropy Probe ( WMAP) and Planck requires a good 
model of the microwave foreground emission to allow the 

[ reliable subtraction of foreground contamination from 
the CMB radiation. 
In the 10-60 GHz frequency range, electric dipole 

I emission (Drainc & Lazarian 1998) from rapidly spin- 
ning tiny dust grains (mostly polycyclic aromatic hy- 
drocarbons, hereafter PAHs) is an important com- 
ponent of Galactic foregrounds that dominates the 
CMB signal (see Planck Collaboration et al. 2013b). 
In the last several years, significant progress has 
been made in understanding spinning dust in terms 
of both theory (Ali-Hai'moud et al. 2009; Hoang et al. 

. 2010; Hoang et al. 2011; Silsbee et al. 2011; sec 
Hoang & Lazarian 2012 for a review) and observation 

■ (Dickinson et al. 2009; Kogut et al. 2011; Tibbs et al. 
2012). In particular, Planck results have con- 
firmed spinning dust emission as the most reli- 
able source of anomalous microwave emission (AME) 
(Planck Collaboration et al. 2011). Planck is poised to 
release interesting results on the CMB polarization; how- 
ever, the question to what extent the spinning dust emis- 
sion contaminates to the polarized CMB signal remains 
unclear. 

The degree of polarization of spinning dust emis- 
sion depends on how efficient the alignment of PAHs 
is in the interstellar magnetic field. Lazarian & Draine 
(2000) (hereafter LDOO) suggested that PAHs can be 
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aligned via resonance paramagnetic relaxation-a mech- 
anism which extends the classical Davis-Greenstein 
(Davis & Grcenstein 1951) mechanism for very fast ro- 
tating grains for which the Barnett magnetization aris- 
ing from fast rotation of grains cannot be neglected. 
They predicted the polarization of spinning dust typi- 
cally < 1% for frequency > 20 GHz. 

Recent observational studies (Lopez-Caraballo et al. 
2011; Maccllari et al. 2011 showed that the upper limit 
for the AME polarization is between 2 — 5%. In addi- 
tion, an upper limit of 1% for the AME polarization is re- 
ported in various media (see Rubiho-Martin et al. 2012). 
Since Planck and other CMB experiments provide ex- 
tremely precise measurements of polarization, a reliable 
prediction for the spinning dust polarization is useful for 
separation of polarized Galactic foreground components 
from the CMB. 

The problem of grain alignment, especially for PAHs'' 
is complicated in general (see Lazarian 2007 for a review), 
but one can derive constraints on grain alignment obser- 
vationally (see Martin 2007). An important attempt to 
obtain observational constraints on grain alignment was 
carried out in a classical paper by Kim & Martin (1995). 
The authors applied maximum entropy method to in- 
fer the mass distribution of aligned grains through fit- 
ting theoretical polarization curves to observation data. 
They discovered that interstellar silicate grains of size 
a > 0.05 /im are aligned, whereas smaller grains, includ- 
ing PAHs, are weakly aligned. Draine & Fraisse (2009) 
performed simultaneous fitting to the typical extinction 
and polarization curve of the diffuse ISM and came to 

^ The accepted theory of grain alignment for large grains 
is that based on radiative torques (Dolginov &; Mitrofanov 
1976; Draine & Weingartner 1996, Draine & Weingartner 
1997; Lazarian & Hoang 2007, Lazarian & Hoang 2008; 
Hoang & Lazarian 2008, Hoang & Lazarian 2009b, 

Hoang & Lazarian 2009a) but the radiative torques are neg- 
ligibly small for PAHs. This encourages us to explore the 
possibilities provided by other alignment mechanisms. 
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the similar conclusion as in Kim & Martin (1995) that 
small (a < 0.05 /xm) grains are weakly aligned and large 
(a > 0.1 fim) grains are efficiently aligned. 

Since the tiny dust grains that radiate spinning dust 
emission are likely the same as those that induce the 
ultraviolet (UV) extinction bump at A = 2175A,'' a nat- 
ural way to search for the alignment of PAHs is through 
its imprint on the UV polarization. While the UV 
extinction bump is ubiquitous in the diffuse interstel- 
lar medium (ISM), the UV polarization bump is rarely 
seen, except for two stars HD 197770 and HD147933- 
4, which show a prominent UV polarization feature 
at A = 2175A (Clayton et al. 1992; Wolff et al. 1993; 
Wolff et al. 1997). In addition to the UV polarization 
excess, HD 197770 exhibits a polarization plateau that 
cannot be accounted for with the typical ISM polariza- 
tion curve-Serkowski law (Serkowski 1973). 

A number of studies have suggested that the UV polar- 
ization bump can arise from aligned small graphite grains 
(Draine 1988; Wolff et al. 1993; Wolff et al. 1997); how- 
ever, a detailed study on the alignment efficiency of small 
graphite grains is not yet available. 

The main goal of this paper is to find the degree of 
alignment of interstellar grains including PAHs that re- 
produces the UV polarization excess seen in HD 197770 
and HD 147933-4 and employ the inferred degree of align- 
ment to predict the degree of polarization of spinning 
dust emission. The paper is structured as follows. 

In §2 we calculate extinction cross section and polar- 
ization cross section for silicate and carbonaceous oblate 
spheroidal grains. In §3, we describe a procedure to de- 
rive grain size distributions and degree of grain alignment 
by fitting theoretical predictions to observed extinction 
and polarization curves and present the obtained results. 
In §4 polarized spinning dust emission is calculated us- 
ing the degree of grain alignment obtained for the best 
fit model. Discussion and summary are presented in §5 
and §6, respectively. 

2. OPTICAL PROPERTIES OF DUST GRAINS 

2.1. Extinction and Polarization Cross Section 

Interstellar dust grains are widely known to induce the 
extinction and polarization of starlight due the absorp- 
tion and scattering of light out of the line of sight by dust 
grains. 

Let us consider an oblate spheroidal grain with the 
symmetry axis ai . A perfectly polarized electromagnetic 
wave with the electric field vector E propagates along 
the z direction, which is perpendicular to the symmetry 
axis. Let Ccxt(E _L ai) and Coxt(E||ai) be the extinc- 
tion of radiation for the cases in which E is parallel and 
perpendicular to ai , respectively. For the sake of simpli- 
fication, we denote these extinction cross sections by C± 
and C||. 

For the general case in which E makes an angle 9 with 
ai, the extinction cross section becomes 

Cext = cos^ 0C\\ + sin^ 0Cj_, (1) 

Since the original starlight is unpolarized, one can com- 
pute the total extinction cross section for a randomly ori- 

Such a UV bump is believed to arise from the electronic transi- 
tion TT — 7r* in the sp^-bonded carbon sheets of small carbonaceous 
grains (see Draine 1989; Draine & Li 2007). 



ented grain by integrating Equation (1) over the isotropic 
distribution of 6*, i.e., fisodO ~ s'm 9 d9. As a result, 

Cext = ^(2Ci+q|). (2) 

The dust extinction cross section can be represented 
through the extinction efficiency, which is defined as the 
ratio of the extinction cross section to the geometrical 
cross section: 

Q,,,(a,A) = ^^^i%^. (3) 

The polarization coefficient for oblate spheroidal grains 
is equal to 

Cpoi - [Cext(E||ai) - Cext(E ± ai)] . (4) 

2.2. Silicate grains 

For silicate grains, we adopt the DDSCAT code 
(Draine & Flatau 2012) to compute the extinction cross 
section C^xt(a, A) using the dielectric functions of amor- 
phous silicate from Draine (2003). 

Figure 1 shows AQext/a and AQpoi/a as functions of 
A^^ for a silicate spheroidal grain with the axis ratio 
a/b ~ 2 and for the different grain sizes. 

2.3. Carbonaceous grains 

Carbonaceous grains consist of a fraction Cpah of PAHs 
and 1 — Cpah of graphite grains with 

CPAH(a) = (1 - ggra) X miu [l, {a>;/af] , (5) 

where = 50A, (/gia ~ 0.01 arc assumed as in 
Draine & Li (2007). 

For graphite grains, which are an anisotropic material, 
wc assume the dielectric function from Draine (2003) 
and compute Cf^t(a, A) using DDSCAT. We consider 
two cases in which E are parallel and perpendicular to 
the grain optical axis (c-axis). Using the ^ — § ap- 
proximation, one can obtain the extinction cross section 
C^^ = (Cext(E||c) + 2Cext(E ± c)) /3 for randomly ori- 
ented graphite grains. 

For PAHs, the extinction C^^^{X) is computed as in 
Equations (5)-(10) in Li & Draine (2001). The polariza- 
tion cross section of PAHs is taken to be the same as 
graphite grains of the same size. 

2.4. Extinction and Polarization Curve 

The extinction induced by randomly oriented grains in 
units of magnitude is defined by 

/ pobs \ 

^(A)-2.51ogi„ ^-j^j = 1.086ta, 

= 1.086 £ Ci,t{a^)nHa,)dz, (6) 

J— carb.sil i—0 

where is the intrinsic flux from the star, F^^^ ~ 
F^e~'^^ is the observed flux, and t\ is the optical depth. 

To find the polarization by aligned grains, let us de- 
fine an observer's coordinate system in which the line of 
sight is directed along the z— axis, and x— and y— axes 
constitute the sky plane. The polarization of starlight 
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Fig. 1. — Extinction efficiency XQcxt/o. (upper) and polarization 
efficiency XQp^i/a (lower) of silicate grains as functions of 1/A for 
the different grain sizes. 

arising from the dichroic extinction by aligned grains in 
a cell of dz is computed as 



dpW = ^ — - 



J2 ^{C.-Cy)n{a,)dz, (7) 

i=0 



where Na is the number of size bin, n{aj) = [dn / da)da is 
the number of grains of size a^, Cx and Cy are the grain 
cross section along x— and y— axes, respectively. 

For the case of perfect internal alignment of grain sym- 
metry axis ai with angular momentum, by transforming 
the grain coordinate system to the observer's coordinate 
system and taking corresponding weights, we obtain 



C,=Ci-%sin2/?, 



(8) 



%i(2cos^/3cos2 7 + sin2/3sin2 7), (9) 



where 7 is the angle between the magnetic field and the 
sky plane and j3 is the angle between the grain angular 
momentum and the magnetic field. 



Cy C'pol ~ 



(3cos2^- 1) 



cos^ 7, 



(10) 



Taking the average of Cx — Cy over the distribution of 
the alignment angle /3, it yields 

Cx - Cy = C^oiQ J cos^l, (11) 

where Qj = {Gj) is the ensemble average of G.j — 
(3cos^/? — 1) /2 that describes the alignment of grain 
angular momentum with the ambient magnetic field. 

When the internal alignment is not perfect, following 
the similar procedure, we obtain 

Cx - Cy = C^oi{GjGx) cos^ 7 = CpoiiJcos^ 7, (12) 

where Gx = (3cos^ 6* — l) /2 with 6 being the angle be- 
tween grain symmetry axis and the angular momentum, 
and R = (GjGx) is the Rayleigh reduction factor (see 
also Roberge & Lazarian 1999). The degree of internal 
alignment is described by Qx = (Gx)- 

For a perpendicular magnetic field, i.e., B lies on the 
sky plane. Equation (12) simply becomes Cx — Cy = 
CpoiR- 

Plugging Equation (12) into this above equation, we 
obtain 



ji— carb.sil i—0 



E i,ClJa,)R^{a,)n\a,)dz, (13) 



where R^ (ui) is the Rayleigh reduction factor for the 
grain specie j of size a^. In the following, the Rayleigh 
reduction factor can also be referred as the alignment 
function. 

Frequently, it is more convenient to represent the ex- 
tinction (polarization) through the extinction (polariza- 
tion) cross section. Hence, the above equations can be 
rewritten as 



A(A)== 1.086 acxt(A) x A^h, 
p(A)==crpoi(A) X Nu, 



(14) 
(15) 



where 7Vh( cm ^) is the column density and (Text and (Tpoi 
in units of cm^H~^ are the dust extinction cross section 
and polarization cross section, respectively. 

3. GRAIN SIZE DISTRIBUTION AND GRAIN ALIGNMENT 
FUNCTION CONSTRAINED BY OBSERVATION 

3.1. Fitting method 

In this section, we find grain size distributions and 
alignment functions by fitting theoretical predictions to 
the observation data for the HD 197770 and HD 147933- 
4 stars. The parameters for these stars, including Eb-v-, 
optical depth Ty, the ratio of visual to selective extinc- 
tion Rv = Av/Eb-v, maximum polarization pmax and 
the wavelength at the peak Amax are listed in Table 1 
(see also Wolff et al. 1997). 

We adopt the mixture dust model consisting of amor- 
phous silicate grains, graphite grains and PAHs (see 
Wcingartner & Draine 2001 ; Draine & Li 2007). We 
note that the grain size distribution may not be unique. 
For example, Draine & Fraisse (2009) obtained four dif- 
ferent models of grain size distribution that can repro- 
duce the diffuse ISM extinction and polarization data. 
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TABLE 1 

Optical extinction and polarization parameters 



Star 


Eb-v 


Rv 


TV 


Amax(Mm) 


Pmax(%) 


HD 197770 
HD 147933-4 


0.58 
0.47 


3.1 
4.3 


1.66 
1.86 


0.511 ±0.001 
0.683 ± 0.003 


3.90 ±0.01 
2.72 ±0.01 



In the present paper, we assume that graphite grains 
of oblate spheroidal shape are randomly oriented. Such 
an assumption allows us to obtain an upper limit 
on the alignment efSciency of PAHs since aligned 
graphite grains would contribute some polarization to 
the 2175A feature. We consider a model in which both 
silicate grains and PAHs are oblate spheroidal. The axis 
ratio a/b — 2 is adopted for the spheroidal grains as in 
the previous section. As in Draine & Fraisse (2009), we 
simultaneously fit theoretical predictions to the extinc- 
tion curve and polarization curve. 

Following Kim & Martin (1995), we find the grain size 
distribution and alignment function by minimizing x^, 
which is given by 



X 



2 I 2 



(16) 



where 



Xoxt = W, 



'ext [^mod(A,) - Aobs(A,)]^ (17) 



1=0 



Xpol = Wpol 



[Pmod(Aj) -pobs(Ai)]^ , (18) 



1=0 



with Wcxt and Wpoi being the fitting weights for the ex- 
tinction and polarization, respectively. Here, the summa- 
tion is performed over N\ wavelength bins. We consider 
Na = 128 size bins from a = 3.6A to 1 iim and Nx = 64 
from A = 0.15 to 2.5 fim. 

In the above equations, the observed extinction Aobs 
is calculated using the formula in Cardelh et al. (1989) 
when Ry is known, and Pobs is obtained from Wolff et al. 
(1997). The extinction Amod(A) and polarization 
J'mod(A) are given by Equations (6) and (13). Follow- 
ing Draine & AUaf-Akbari (2006), we also introduce the 
constraints for the non-smoothness of the grain size dis- 
tribution and alignment function. 

The fitting procedure is started with an initial size 
distribution n(a) that best reproduces the observational 
data for the diffuse ISM, which corresponds to model 3 
in Draine & Fraisse (2009). By doing so, we implicitly 
assume that dust properties are similar throughout the 
ISM and the difference in the polarization of starlight is 
mainly due to the efficiency of grain alignment, which de- 
pends on environment conditions along the line of sight, 
e.g., radiation field, magnetic fields and gas density. For 
silicate grains, we take the alignment function for the 
ISM from Draine & Fraisse (2009) as an initial alignment 
function. For PAHs, we consider two cases. In the first 
case, we take the initial alignment function to be equal to 
that of silicate grains. In the second case, we assume an 
initial alignment function with a peak at a = lOA from 
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Fig. 2. — Observed extinction curve (filled ellipses) versus model 
for HD 197770 and HD 147933-4. An excellent fit achieved for HD 
197770. 

the theoretical prediction based on resonance paramag- 
netic alignment (see Hoang et al. 2013). 

The fitting is carried out by the least square method 
based on Monte Carlo simulations. Basically, for each 
size bin, we add a small amount 6R and Sn generated 
from a uniform distribution to R{a) and n{a). Then, we 
calculate Amod and Pmod for the new degree of alignment 
and size distribution using Equations (6) and (13). The 
value of is also calculated. This sampling process is 
iterated for A'rand = 100 times with each size bin to find 



the 



minimum x 



3.2. Results 



Figure 2 shows the extinction curves from our model 
compared to the observed data for the HD 197770 and 
HD 147933-4 stars. The model can reproduce the ob- 
served extinction of HD 197770 very well. For HD 
147933-4, the model underestimates the extinction for 
A = 0.17- 0.2 /.tm. 

Figure 3 shows the mass distributions of silicate and 
carbonaceous grains for the best fit model for HD 197770 
(upper) and HD 1479333-4 (lower). The distribution 
functions (DL07 and DF09) from Draine & Li (2007) and 
Draine & Fraisse (2009) are shown for comparison. 

For HD 197770 with Ry = 3.1, it can be seen that 
the mass of PAHs in our model is the same as DL07 
and DF09. The mass distribution of graphite grains is 
slightly different from DF09 for a = 0.04 - 0.08 ^m, but 
it is distinct from DL07. For silicate grains, the mass 
distribution peaks at a w 0.07 /im. Compared to DF09, 
our model requires a slight increase of silicate mass to 
account for the polarization plateau seen in HD 197770. 

For HD 147933-4 with Ry = 4.3, the mass distribution 
of silicate grains peaks at a larger size (a w 0.25 fim) than 
for HD 197770, whereas the mass of 0.03 — 0.2 fim silicate 
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Fig. 3. — Mass distribution of silicate grains and carbonaceous 
grains for the best fit model (solid lines) for the HD f 97770 (upper) 
and HD f47933-4 (lower) stars. Dotted (DF09) and dashed (DL07) 
lines show the grain mass distributions from Drainc & Fraisse 
(2009) and Draine & Li (2007). 

grains is reduced. The peak in the mass of 0.08 fim) 
graphite grains is also increased compared to DF09. Such 
a decrease in the mass of small grains and increase in the 
mass of large grains is necessary to reproduce the value 
of Rv higher than the typical value {Ry = 3.1) of the 
diffuse ISM. 

Figure 4 shows the best fit model versus the observed 
polarization data from Wolff et al. (1997). Here we show 
the error bars of 3cr. For HD 197770, it can be seen that 
the model reproduces the 2175 A polarization feature 
very well. The UV polarization bump in HD 147933-4 
can also be accounted for, but the model shows a broad 
bump from = 4.0 — 5.8 fim~^ around the 2175A fea- 
ture. 

Figure 5 (upper) shows the alignment function R{a) 
that reproduces the polarization curve in Figure 4 for HD 
197770. We can see that the alignment of silicate grains 
rises rapidly from a = 0.06 fim and becomes constant 
with R = 0.46 for a > 0.2 /im. For a < 0.06 /im, the de- 
gree of alignment declines more or less smoothly with the 
increasing a but is still considerable (i?(a) = 0.01 — 0.05) 
for very small silicate grains of a = 0.003 — 0.06 ^m. The 
alignment of a < 0.003 /xm (i.e. 30A) silicate grains is 
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Fig. 4. — Polarization curves for our model versus the observed 
polarization data for HD f 97770 (upper) and HD 147933-4 (lower). 
The 3a- error bars are shown. 

neglected because the polarization by tiny silicate grains 
containing minor dust mass is negligible. 

In particular, one can see from the figure that the best 
fit model corresponds to a peaky alignment function of 
PAHs with i?peak ~ 0.014. Moreover, a rather low align- 
ment degree of PAHs could be sufficient to reproduce 
the UV polarization feature since the mass of PAHs is 
dominant for the dust mass for a < lOA ( sec Figure 3). 

For HD 147933-4 (see Figure 5, lower), we find that the 
best fit requires no alignment of PAHs (i.e., the align- 
ment of only silicate grains can reproduce the observed 
polarization). The degree of alignment R = 0.72 for 
a > 0.15 /im silicate grains. Moreover, the alignment of 
silicate grains decline steeply for a < 0.1 and becomes 
negligible (i? < 0.01) for a < 0.05 /im. It turns out that 
small silicate grains a < 0.05 fira have minor contribution 
to the polarization, which is different from HD 197770. 
The sharp decline in alignment of small (a < 0.1 ^m) sil- 
icate grains (i.e., the polarization is dominated by larger 
silicate grains) is required to account for the large value 
of the peak wavelength Amax seen in HD 147933-4. 

4. POLARIZATION OF SPINNING DUST EMISSION 
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Fig. 5. — Upper panel: Degree of alignment as a function of 
grain size for silicate grains (blue line) and carbonaceous grains 
(red line) obtained from the fitting for HD 197770. Dotted line 
shows a = 3 X 10~^cm from which the polarization by aligned 
small silicate grains becomes negligible. Lower panel: Similar to 
upper panel but for HD 147933-4. The alignment of silicate grains 
is shown only due to the lack of PAH alignment (see the text). 

4.1. Degree of grain alignment for spinning dust 
emission 

To calculate polarized spinning dust emissivity, it is im- 
portant to understand the relation between the degree of 
alignment of PAHs that is responsible for the polariza- 
tion in microwave emission and the degree of alignment 
of PAHs that results in starlight polarization by extinc- 
tion (i.e., Rayleigh reduction factor). 

As discussed in Lazarian (2007), these two measures 
are expected to be different. Indeed, even spherical PAHs 
aligned by paramagnetic mechanism can produce polar- 
ized microwave emission. Naturally, UV absorption or in- 
frared emission is not expected from the spherical PAHs. 
Although we do not really believe that most of PAHs are 
spherical, this illustrates the problem that we face com- 
paring different observational consequences of alignment. 
For instance, fast wobbling and and flipping grains are 
also expected to deliver substantially reduced degrees of 
polarization in terms of UV absorption. However, this 
just decreases the polarization of microwave emission by 
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Fig. 6. — Alignment measures as functions of grain size a: 
Rayleigh reduction factor R, degree of alignment of grain axes with 
the angular momentum Qxi and degree of alignment of angular 
momentum with the magnetic field Qj. 

a factor of unity''. 

In the following, to find the upper limit of polarized 
spinning dust, we adopt the degree of alignment of angu- 
lar momentum with the magnetic field Qj. Thus, from 
the inferred alignment function R(a), we can derive Qj 
using the following relationship (see Section 2.4) 

R{a) = (GxGj) = Qj X Qx{l + /corr), (19) 

where Qx describes the degree of alignment of grain 
axes with the angular momentum and /corr describes the 
correlation between the internal and external alignment. 
The case /corr = indicates that the internal alignment 
is independent from the external alignment. Our calcu- 
lations shows that /corr ^ 0.6 for the range of size of 
PAHs. 

Figure 6 shows the alignment measures as functions 
of grain size. The Rayleigh reduction factor R is from 
the best fit model (Figure 5, upper). Qx is obtained 
from detailed calculations in Hoang et al. (2013) using 
Langevin equations for the cold neutral medium (CNM) 
and a constant dust temperature = 60 K. Qj is cal- 
culated using Equation (19). It can be seen that Qj is a 
factor of 3 higher than the Rayleigh reduction factor. 

4.2. Polarized spinning dust emissivity 

The polarized emissivity and unpolarized emissivity of 
spinning dust emission can be given by 



dn 



1'^ = J Qj(a).7iy(a)— da, 
Ji. = y".7i.(a)^da, 



(20) 
(21) 



where (a) is the spinning dust emissivity at frequency 
1/ from a grain of size a. The frequency dependence of po- 
larization of spinning dust emission is p{h') = <lvl jv We 

^ One can say that for the microwave polarization from PAHs 
the alignment in terms of the PAH angular momenta is impor- 
tant, while for the UV absorption polarization, it is the alignment 
in terms of grain axes that is essential. The two alignments are 
different (see Lazarian & Robergc 1997;Hoang & Lazarian 2008; 
Hoang et al. 2011). 
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V (GHz) 

Fig. 7. — Upper panel: Unpolarized (solid line) and polarized 
(dotted line) spinning dust cmissivity from the diffuse ISM along 
the line of sight to HD 197770. Lower panel: Frequency dependence 
of degree of polarization of spinning dust emission. 

adopt the spinning dust model from Hoang et al. (2011) 
for oblate spheroidal grains with the ratio axis ajh — 2. 

Figure 7 (upper panel) shows the total (unpolarized) 
and polarized spinning dust emissivity for the best fit 
model. The degree of polarization is shown in Figure 7 
(lower panel). It can be seen that the maximum degree 
of polarization is about 3.2% and peaks at « 5 GHz. 
The polarization is < 1.8% for v > 20 GHz. 

5. DISCUSSION 

In the Planck era, electric dipole emission from rapidly 
rotating PAHs has become an accepted component of dif- 
fuse galactic foreground emissions that contaminate to 
the CMB radiation. In light of more exciting Planck re- 
sults on polarization, a pressing issue remains is to quan- 
tify the level of contamination by polarized spinning dust 
emission to the polarized CMB signal. The present pa- 
per seeks constraints on the polarization of spinning dust 
emission based on the models that best fit to observed 
polarization data of the stars having potential evidence 
of PAHs alignment. 

5.1. UV polarization excess and alignment of PAHs 



The polarization excess at 2175A toward two stars, 
HD197770 and HD147933-4, was discovered a long time 
ago (Clayton et al. 1992; Wolff et al. 1993), and it is 
widely believed that the excess originates from aligned 
PAHs. Employing the observation data for these stars, 
we obtain the grain size distributions and alignment func- 
tions that best reproduce both the extinction and polar- 
ization curves. 

For HD 197770, we found that the alignment of sili- 
cate grains only cannot reproduce the UV polarization 
bump adequately. Instead, a model with aligned silicate 
grains plus weakly aligned PAHs can successfully repro- 
duce the UV polarization bump as well as the polariza- 
tion plateau. The inferred degree of PAHs alignment 
present in HD 197770 varies with the grain size and has 
peak of 0.014 at a « 9A. Ahhough the degree of PAH 
alignment is rather small, due to the dominance of PAHs 
for the dust mass for a < 20A, it is sufficient to repro- 
duce the UV polarization excess. In fact, the low degree 
of alignment for PAHs is not unexpected from theoretical 
predictions based on resonance paramagnetic alignment, 
which was proposed by Lazarian & Draine (2000) and 
numerically studied in Hoang et al. 2013. 

For HD 147933-4, on the other hand, we found no in- 
dication of alignment for PAHs. Yet the model with only 
aligned silicate grains can account for the UV polariza- 
tion excess; although the model predicts a broad bump 
(see Figure 5).^ 

Before understanding why there is such a difference 
in two stars, let us recall the significant difference in 
optical and polarization properties of HD 197770 and 
HD 147933-4 (see Table 1). Indeed, the latter has a 
much higher ratio of visual-to-selective extinction Ry 
and a much larger peak wavelength Amax- The polariza- 
tion ratio p(2175A)/p„iax is much lower in the later case. 
Qualitatively, Ry reflects the average size of dust grains, 
whereas Amax reflects the average size of aligned grains. 
Thus, compared to HD 197770, HD 147933-4 essentially 
has the larger average grain size and the larger average 
size of aligned grains. As a result, the UV polarization is 
dominated by the alignment of large grains. Moreover, 
since the large silicate grains have optical properties with 
oscillating features in the UV due to interference effects,^ 
they can induce the broad polarization feature seen be- 
yond A~^ « 4.6 ^m~^. 

5.2. Why is the UV polarization excess not seen for 
most stars? 

The question now is that if the PAHs are potentially 
aligned, why we do not see the UV polarization feature 
for most of stars. That answer seems straightforward. 
Indeed, let us consider the lines of sight of similar Ry 
and Amax- Then, the presence of UV polarization excess 
depends on the degree of alignment of PAHs in the ambi- 
ent magnetic field. Although the alignment mechanism 
for PAHs is still not clear, the most promising mecha- 
nism is based on the resonance paramagnetic relaxation 
(LDOO). Since the paramagnetic alignment is sensitive to 

^ More discussions on why it is very challenging to obtain a 
satisfactory fit for the UV polarization data for the stars with 
Amax > 0.6 /^m can be found in Kim & Martin (1995). 

^ The oscillating feature is prominent for dust with the dielectric 
function having small imaginary part (weakly absorbing material) . 
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the magnetic field and gas randomization, the variation 
of the magnetic field (both strength and direction) and 
gas density along different lines of sight can result in the 
absence of the UV polarization feature. 

In addition, the resonance paramagnetic alignment de- 
pends on spin-spin and spin-lattice relaxation within 
PAHs. If the relaxation is reduced, the expected de- 
gree of grain alignment is also reduced. The point of the 
relaxation in PAHs, as it was claimed in LDOO can be 
settled via laboratory studies. 

Finally, the alignment degree of PAHs in the magnetic 
field in general is rather weak (i? < 0.014) for which 
the resulting UV polarization excess would be small. As 
pointed out in Wolff et al. (1997), this feature might not 
be detected by low signal-to-noise observations. 



5.3. Constraint on the polarization of spinning dust 
emission 

We calculated the degree of polarization for spinning 
dust emission using the alignment function for the best 
fit model to the HD 197770 data. We found that the 
upper limit for the polarization of spinning dust emission 
is ~ 3.2% at = 5 GHz and the degree of polarization 
declines rapidly to below 1.8% for i' > 20 GHz. 

It is worth noting that our calculations assume that the 
UV polarization bump is produced by PAHs only (i.e., 
the contribution of large graphite grains is disregarded). 
In fact, if graphite grains can be aligned, then their con- 
tribution to the UV polarization bump would reduce our 
estimate for spinning dust polarization. Therefore, this 
constraint is indeed the upper constraint for the spinning 
dust polarization. 

Observationally, a number of studies indicate that 
anomalous microwave emission (AME) is weakly polar- 
ized (Mason et al. 2009; Lopez-Caraballo et al. 2011). 
Using WMAP data, Macellari et al. (2011) found the up- 
per limit for the polarization of AME from 1.4 — 2% for 

= 20 - 40 GHz. Rubifio-Martfn et al. (2012) reviewed 
in great detail on observational constraints for the AME 
polarization and an upper limit of 1% for = 20—30 GHz 
is observed for the various environments. It appears that 
our upper limit on spinning dust polarization is consis- 
tent with the current observation data for the polariza- 
tion of AME. 



degree of dust polarization as follows: 

^A.poi ^ q'poi(A) 
-^A a'cxt(A)' 



P., 



(22) 



5.4. Implication for polarized far-infrared thermal dust 

emission 

Planck Collaboration et al. (2013a) shows that the de- 
gree of polarization of thermal dust emission can reach a 
high level Pom > 15% ai v — 353 GHz. Let us estimate 
the degree of dust polarization at this frequency along 
the lines of sight toward two stars HD 197770 and HD 
147933-4. 

For simplicity, we assume that all grains have the same 
temperature. Thus, using the grain size distributions and 
alignment functions from the model, we can obtain the 



where /a,po1 is the polarized emission, I\ is the total dust 
emission, and (Text and cTpoi are present in Equations (14) 
and (15). 

At = 353 GHz (i.e., A ~ 850 //m), wc found the 
degree of polarization Pom = 9.9% and 14.5% for HD 
197770 and 147933-4, respectively. The higher degree 
of polarization for the HD 147933-4 star is expected be- 
cause the degree of alignment of large (a > 0.15 ^m) 
grains which dominate far-infrared dust emission is a fac- 
tor of 1.55 higher than that for the HD 197770 star (see 
Figure 5). It can be seen that a high level of polariza- 
tion (> 15%) seen by Planck is expected from our model 
(- 14.5%). 

To study a potential correlation between the degree 
of polarization of far-infrared dust emission to that by 
dust extinction, we introduce a correlation parameter 
r, which is equal to the ratio of Pom at v = 353 GHz 
to Pcxt/T at A = Ainax- We obtain r « 4.7 for HD 
197770 and r « 7.5 for HD 147933-4. Interestingly, the 
difference in the correlation factor r of the two stars, 
r(HD197770)/r(HD147933 - 4) = 1.59, is close to the 
ratio of the degree of alignment for large grains obtained 
for the two stars. 

6. SUMMARY 

In the present paper, we obtain the constraint on the 
degree of polarization of spinning dust emission based 
on the observed extinction and polarization curves of 
starlight. Our principal results are summarized as fol- 
lows: 

1. A model of grain size distribution and alignment 
function for interstellar grains that best fits to the 
observed extinction and polarization curve is ob- 
tained for two stars HD 197770 and HD147933-4 
with the prominent polarization features at 2175A. 

2. We find that a small degree of alignment (« 0.014) 
of PAHs included to the aligned silicate grains is 
sufficient to reproduce the UV polarization excess 
for HD 197770. For HD 147933-4, we find no indi- 
cation of alignment for PAHs and the alignment of 
silicate grains can account for the excess. 

3. We calculate the polarized spinning dust emis- 
sion using the alignment function from the best fit 
model. We show that the degree of polarization 
for spinning dust emission has a peak of 3.2% at 
ly = 5 GHz and rapidly decline to below 1.8% for 
J/ > 20 GHz. 

4. The degree of polarization for thermal dust emis- 
sion at 353 GHz is estimated to be Pem = 9.9% for 
the HD 197770 star and Pen ^ 14.5% for the HD 
147933-4 star. 
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